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ABSTRACT 

We present an absorption-line survey of optically thick gas clouds - Lyman Limit Systems (LLSs) 

- observed at high dispersion with spectrometers on the Keck and Magellan telescopes. We measure 
column densities of neutral hydrogen Vhi and associated metal-line transitions for 157 LLSs at 2 ;lls = 

1.76 — 4.39 restricted to cm“^ < Vhi < 10^° '^cm“^. An empirical analysis of ionic ratios 

indicates an increasing ionization state of the gas with decreasing and that the majority of LLSs 
are highly ionized, confirming previous expectations. The Si"'"/H*’ ratio spans nearly four orders-of- 
magnitude, implying a large dispersion in the gas metallicity. Fewer than 5% of these LLSs have no 
positive detection of a metal transition; by 2 ~ 3, nearly all gas that is dense enough to exhibit a very 
high Lyman limit opacity has previously been polluted by heavy elements. We add new measurements 
to the small subset of LLS (« 5 — 10%) that may have super-solar abundances. High Si+/Fe''" ratios 
suggest an a-enhanced medium whereas the Si“'"/C"'" ratios do not exhibit the super-solar enhancement 
inferred previously for the Lya forest. 

Subject headings: absorption lines - intergalactic medium — Lyman limit systems 


1. INTRODUCTION 

As a packet of ionizing radiation (hv > I Ryd) tra¬ 
verses the universe, it has a high probability of encoun¬ 
tering a slab of optically thick, H I gas. For sources in 
the z ^ 4 universe the mean fr ee path is only « 30 Mpc 
(physical; IWorseck et al.||2014D, i.e. less than 2% of the 
event horizon. (Jbservationally, researchers refer to this 
optically thick gas as Lyman limit systems (LLSs) owing 
to their unmistakable signature of continuum opacity at 
the Lyman limit (« 9I2A) in the system restframe. A 
fraction of this gas lies within the dense, neutral interstel¬ 
lar medium (ISM) of galaxies, yet the majori ty of opac- 


Lombardi et al.||l994l 

samples. 


For many decades , LLS have been surveyed in quasar 
spectra (e.g. Tytl^|1982 Sargent et al. 1989 Storrie- 


albeit often from heterogeneous 


These works established the high incidence 
of LLSs which evolves rapidly with redshift. With the 
realization of massive spectral datasets, a renaissance 
of LLS surveys has followed yielding statistically ro¬ 
bust measurem ents from homogenous and well-selected 
quasar samples (Prochaska et al.)20I0 Songaila & Cowie 
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20101 IRibaudo et al.||20111 |0’Meara et al. 

[2M3l 

Fuma-| 

galli et al. 

2013 

. Analysis of these hundreds of systems 
ence of approximately 1.2 systems per 

reveals an 

incid 


unit redshift at z ~ 3 that e volves steeply with redshift 


z) PC (iT-z)*'^ for z « 1 — 5 (Ribaudo et al. 


2011 


Fuma- 


galli et aL]|20I3). With these same spectra, researchers 


have further r neasured the mean free path of ionizing ra¬ 
diation (A^fp; [Prochaska et al.|2009 O’Meara et al. 2013 


et al.||2011b| Ribaudo et al. 

2011 

). Indeed, the interplay 


between galaxies and the L 
research which includes stu 

Ly is a highly active area of 
dies of the so-called circum- 

galactic medium (CGM; e.j 

Steidel et al. 

2010 Werk 


et al.l2013 

Prochaska et al. 

(2014a|. - 



Fumagalli et al.|2013||Worseck et al.|2014p, which sets the 

intensity and shape of the extragalactic UV background. 
Following the redshift evolution of the LLS incidence, 
also evolves steeply with the expanding universe, 
implying a more highly ionized u niverse with advancing 
cosmic time (Worseck et al. 2014). 


The preponderance of LLSs bespeaks a major reservoir 
of baryons. In particular, given the apparent paucity of 
heavy elements wit hin galaxies (e.g. Bouche et al.||200^ 


Peeples et al.|2014 ), the LLSs may present the dorninant 


reserv oir ot nietals in the universe (e.g. Prochaska et al.j 
2006|. However, a precise calculation ot the heavy ele- 


ments within LLSs and their contribution to the cosmic 
budget has not yet been achieved. Despite our success at 
surveying hundreds of LLSs, there have been few studies 
resolving their physical properties and thes e have gener¬ 


ally examined a few individual cases ( e.g. Steidel 1990 


Proch aska 1999) or composite spectra (Fumagalli et al. 


2013). 'i'ms reflects both the challenges related to data 


acquisition and analysis together with a historical focus 
in the community towards the ISM of galaxies (probed by 
DLAs) and the more diffuse intergalactic medium (IGM). 

At z > 2, a few works have examined the set of LLSs 
with high H 1 column density (A^hi > 10*^®cm“^), gener¬ 
ally termed the super-LLSs or sub-damped Lya systems. 
Their A^hi frequency distribution /(A^hiW) chemi- 


sample (JDessauges-Zavadsky et al.||2003 |Peroux et al. 

2005 

lO’Meara et al.||20071 |Zafar et al.| 

20131 Som et al. 

2013 

). Ignoring ionization corrections 

whicn may not 
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be justified, these SLLSs exhibit metallicities of approx¬ 
imately 1/10 solar, comparable to the enrichm ent level 
of the highe r- -/Vht, damped Lya systems (DLAs; Rafelski 
et al.||2012 ). In addition, a few LLSs have received spe¬ 
cial attention owing to their peculiar metal-enr ichment 


(Prochaska et al. 


2006 Fumagalli et al. 2011a) and/or 


the detecti on of D tor ^dies ot Pig Bang JNucleosyn - 


thesis (e.g. |Burles & Tytler||1998| fO’Meara et al. 2006|). 
Most recently, a sample of 15 lLs s has been surveye d 
for highly ionized O VI absorption (Lehner et al. 2014), 
which is present at a high rate. A comprehensive study 
of the absorption-line properties of the LLSs at high red- 
shift, however, has not yet been performed. 

Scientifically, we have two primary motivations to sur¬ 
vey the LLSs at z > 2. First and foremost, we aim to 
dissect the physical nature of the gas that dominates the 
opacity to ionizing radiation in the universe. One sus¬ 
pects that these LLSs trace a diverse set of overdense 
structures ranging from galactic gas to the densest fila¬ 
ments of the cosmic web. Such diversity may manifest 
in an wide distribution of observed properties (e.g. metal 
enrichment, ionization state, kinematics). Second, mod¬ 
ern theories of galaxy formation predict that the gas fu¬ 
eling star for mation accretes on t o galaxies in cool, dense 


streams (e.g. Keres et al. 


diative transfer analysis o 


2005 Dekel et al. 2009). Ra- 
hydrodynamic simulations of 


this process predict a relatively high cr oss-section of op- 


&: KeresI 120111 IFumagalli et al.||2011b 

\Mi\ [J^'aucW- 

Giguere et ai.|2014|). Ino 

leed, an optically thick CGM en- 
cies hosting z ~ 2 quasars (Hen- 

velops the massive gala? 

nawi & Prochaska|2007 

Prochaska et al. 12013), LLSs are 


and such gas persists around present -day L"" galaxies (e.g. 


Chen et al. 2010 Werk et al. 2013). The latter has in- 
spired, m part, surveys of the LLbs at z < 1 with ul¬ 


traviolet sp ectroscopy (e.g. |Ribaudo et al.||2011t |Lehner 

eTaLjlM^ . 


THus motivated, we have obtained a large dataset of 
high-dispersion spectroscopy on z > 3 quasars at the 
Keck and Las Campanas Observatories. We have sup¬ 
plemented this program with additional spectra obtained 


to survey the damped Lya systems (e.g. Prochaska et al. 


2007 Berg et al.|2014|) and the intergalactic medium (e.g. 


Faucher-Giguere et al. 2008). In this paper, we present 
the comprehensive dataset of column density measure¬ 
ments on over 150 LLSs. Future manuscripts will ex¬ 
amine the metallicity, chemical abundances, kinematics, 
and ionization state of this gas. This manuscript is out¬ 
lined as follows: Section [^describes the dataset analyzed 
including a summary of the observations and procedures 
for generating calibrated spectra. We define an LLS in 
Section and detail the procedures followed to estimate 
the H I column densities in Section Section presents 
measurements of the ionic column densities am the pri¬ 
mary results of an empirical assessment of these data are 
given in Section A summary in Section concludes 
the paper. 

2. DATA 

This section describes the steps taken to generate a 
large dataset of high-dispersion, calibrated spectra of 
high redshift LLSs. 


2.1. Our Survey 

The sample presented in this manuscript is intended 
to be a nearly, all-inclusive set of LLSs discovered in the 
high-dispersion (echelle or echellette; R > 5,000) spec¬ 
tra that we have gathered at the Keck and Magellan 
telescopes. Regarding Keck, we have examined all of 
the data obtained by Principal Investigators (Pis) A.M. 
Wolfe and J.X. Prochaska at the W.M. Keck Observa¬ 
tory through April 2012, and from Pis Buries, O’Meara, 
Bernstein, and Fumagalli at Magellan through July 2012. 
We also include the Keck spectra analyzed by IPenprasel 

erall(l2^. 

Each spectrum was visually inspected for the presence 
of damped Lya absorption and/or a continuum break at 
wavelengths A < 912A in the quasar rest-frame. The 
complex combination of spectral S/N, wavelength cover¬ 
age, and quasar emission redshift Zem leads to a varying 
sensitivity to an LLS. No attempt is made here to define 
a statistical sample, e.g. to assess the random incidence 
of LLSs nor their A^hi frequency distribution /(TVhd A)- 
We refer the reader to previous manuscripts on this topic 


(Prochaska et al. 2010 Fumagalli et al. 


2013). Because 


our selection is based solely on H I absorption, however 
we believe the sample is largely unbiased with respect 
to other properties of the gas, e.g. metal-line absorption, 
kinematics, ionization state. 

The sample was limited during the survey by: (1) 
generally ignoring LLSs with absorption redshifts within 
3000 km s“^ of the reported quasar redshift Zem, so-called 
proximate LLS or PLLS; and (2) generally ignoring LLSs 
with A^hi < lO^^-^cm”^, especially when the S/N was 
poor near the Lyman limit. We note further that many 
of the Keck spectra were obtaine d to study damped Lya 


systems (DLAs) at z > 2 (e.g. Prochaska et al. 


20071 [Rafelski et al.||2012[ |Neeleman et al.||2013 


2001 

BergI 


et al.||2014 ). We have ignored systems targeted as f)LAs 
and also absorbers within « 1500km s“^ of these DLAs 
because the DLA system complicates analysis of the H I 
Lyman series and metal-line transitions of any nearby 
LLS. In § [^ we offer a strict definition for an LLS to 
define our sample of 157 systems. 

2.2. Observations 

We present data obtained at the W.M. Keck and Las 
Campanas Observatories using the twin 10 m Keck I and 
Keck II telescopes and the twin 6.5 m Baade and Clay 
telescopes. Altogether, we used four spectrometer s: (1) 


et al. 

1994): (2) the Echellette Spectrograph and Imager 

(ESI; 

Sheinis et al.[[2002): (3) the Magellan Inamori Ky- 

ocera 

Echelle (MIKE; 

Bernstein et al. 12003): and (4) the 

iVlagellan Echellette Spectrograph (IVlagE; 

Marshall et al. 


2008). 


The MagE spectra were presented in [Fumagalli et al.j 


(2013) and we refer the reader to that manuscript for 
details on the observations and data reduction. Similarly 
the ESI observat ions have been published previously in a 


series of pai 

pers ([Prochaska et al. [[2003b 

2007 

O’Meara 

et al.[[2007 

Ratelski et al.||2U12p. 


Observing logs for the 
provided in Tablesand [^ A significant f raction of these 
data have been anmyzedpreviously (e.g. O’Meara et al. 


2007 Eaucher-Giguere et al.|2008 Neeleman et al.|2013 ( 
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TABLE 1 

JOURNAL OF HIRES OBSERVATIONS 


QSO 

Alt. Name 

RA 

(J2000) 

DEC 

(J2000) 

r/V" 

(mag) 

^em 

Date 

(UT) 

Slit'’ 

Mode 

Exp 

S/N'’ 

(pix-i) 

SDSS0121+1448 


01;21:56.03 

+14:48:23.8 

17.1 

2.87 

08 Sep 2004 

Cl 

HIRESb 

7200 

15/26 

PSS0133+0400 


01:33:40.4 

+04:00:59 

18.3 

4.13 

27 Dec 2006 

Cl 

HIRESr 

7200 

14/20 

SDSS0157-0106 


01:57:41.56 

-01:06:29.6 

18.2 

3.564 

18 Dec 2003 

C5 

HIRESr 

9000 

X/U 

Q0201+36 


02:04:55.60 

+36:49:18.0 

17.5 

2.912 

06 Oct 2004 

Cl 

HIRESb 

3600 

4.5/9 

PSS0209+0517 


02:09:44.7 

+05:17:14 

17.8 

4.18 

18 Sep 2007 

Cl 

HIRESr 

11100 

31/24 

Q0207-003 


02:09:51.1 

-00:05:13 

17.1 

2.86 

08 Sep 2004 

Cl 

HIRESb 

5400 

15/40 







09 Sep 2004 

Cl 

HIRESb 

8100 


LB0256-0000 


02:59:05.6 

+00:11:22 

17.7 

3.37 

03 Jan 2006 

C5 

HIRESb 

7049 

11/17 

Q0301-005 


03:03:41.0 

-00:23:22 

17.6 

3.23 

09 Sep 2004 

Cl 

HIRESb 

7800 

Xh5 

Q0336-01 


03:39:01.0 

-01:33:18 

18.2 

3.20 

26 Oct 2005 

C5 

HIRESb 

3600 

X/10 







01 Nov 2003 

Cl 

HIRESr 

10800 

15 

SDSS0340-0159 


03:40:24.57 

-05:19:09.2 

17.95 

2.34 

06 Oct 2008 

Cl 

HIRESb 

3000 

7/15 

HE0340-2612 


03:42:27.8 

-26:02:43 

17.4 

3.14 

26 Oct 2005 

Cl 

HIRESb 

7200 

17/X 

SDSS0731+2854 


07:31:49.5 

+28:54:48.6 

18.5 

3.676 

04 Jan 2006 

C5 

HIRESb 

7200 

X/15 

Q0731+65 


07:36:21.1 

+65:13:12 

18.5 

3.03 

28 Oct 2005 

C5 

HIRESb 

5400 

X/16 







04 Jan 2006 

C5 

HIRESb 

7200 

X/12 

J0753+4231 


07:53:03.3 

+42:31:30 

17.92 

3.59 

26 Oct 2005 

C5 

HIRESb 

3300 

X/12 







28 Oct 2005 

C5 

HIRESb 

4800 

X/16 

SDSS0826+3148 


08:26:19.7 

+31:48:48 

17.76 

3.093 

27 Dec 2006 

Cl 

HIRESr 

7900 

ZllTl 

J0828+0858 


08:28:49.2 

+08:58:55 

18.30 

2.271 

14 Apr 2012 

Cl 

HIRESb 

1295 

6/9 

J0900+4215 


09:00:33.5 

+42:15:46 

16.98 

3.290 

15 Apr 2005 

Cl 

HIRESb 

4700 

X/20 

J0927+5621 


09:27:05.9 

+56:21:14 

18.22 

2.28 

14 Apr 2005 

C5 

HIRESb 

8500 

6/20 

J0942+0422 


09:42:02.0 

+04:22:44 

17.18 

3.28 

18 Mar 2005 

Cl 

HIRESb 

7200 

27/X 

J0953+5230 


09:53:09.0 

+52:30:30 

17.66 

1.88 

18 Mar 2005 

Cl 

HIRESb 

7200 

18/22 

Q0956+122 


09:58:52.2 

+12:02:44 

17.6 

3.29 

03 Jan 2006 

C5 

HIRESb 

7200 

X^O 







07 Apr 2006 

Cl 

HIRESr 

1800 

15/10 

HSlOll+4315 


10:14:47.1 

+43:00:31 

16.1 

3.1 

14 Apr 2005 

C5 

HIRESb 

5100 

X/40 







27 Apr 2007 

B2 

HIRESr 

3600 

47/47 







28 Apr 2007 

B2 

HIRESr 

3600 


J1019+5246 


10:19:39.1 

+52:46:28 

17.92 

2.170 

11 Apr 2007 

Cl 

HIRESb 

7200 

11/16 

Q10I7+109 


10:20:10.0 

+10:40:02 

17.5 

3.15 

06 Apr 2006 

C5 

HIRESb 

7200 

25/X 

J1035+5440 


10:35:14.2 

+54:40:40 

18.21 

2.988 

25 Mar 2008 

Cl 

HIERSr 

10800 

23/24 

SDSS1040+5724 


10:40:18.5 

+57:24:48 

18.30 

3.409 

04 Jan 2006 

C5 

HIRESb 

8100 

xh2 

Q1108-0747 


11:11:13.6 

-08:04:02 

18.1 

3.92 

07 Apr 2006 

Cl 

HIRESr 

7200 

30/10 

J1131+6044 


11:31:30.4 

+60:44:21 

17.73 

2.921 

26 Dec 2006 

Cl 

HIRESb 

7200 

14/18 

J1134+5742 


11:34:19.0 

+57:42:05 

18.20 

3.522 

05 Jan 2006 

C5 

HIRESr 

6300 

26/22 

J1159-0032 


11:59:40.7 

-00:32:03 

18.10 

2.034 

14 Apr 2012 

Cl 

HIRESb 

2400 

5/7 

Q1206+1155 


12:09:18.0 

+09:54:27 

17.6 

3.11 

06 Apr 2006 

C5 

HIRESb 

7200 

23/X 

Q1330+0108 


13:32:54.4 

+00:52:51 

18.2 

3.51 

07 Apr 2006 

Cl 

HIRESr 

7200 

11/9 

HS1345+2832 


13:48:11.7 

+28:18:02 

16.8 

2.97 

14 Apr 2005 

C5 

HIRESb 

4800 

X/27 

PKS1354-17 


13:57:06.07 

-17:44:01.9 

18.5 

3.15 

28 Apr 2007 

C5 

HIRESr 

7200 

%ll 

J1407+6454 


14:07:47.2 

+64:54:19 

17.24 

3.11 

14 Apr 2005 

C5 

HIRESb 

5400 

X/20 

HS1431+3144 


14:33:16.0 

+31:31:26 

17.1 

2.94 

06 Apr 2006 

C5 

HIRESb 

6000 

25/43 

J1454+5114 


14:54:08.9 

+51:14:44 

17.59 

3.644 

14 Jul 2005 

C5 

HIRESr 

1800 

10/7.5 

J1509+1113 


15:09:32.1 

+ 11:13:14 

19.0 

2.11 

15 Apr 2012 

Cl 

HIRESb 

5200 

4/7 

J1555+4800 


15:55:56.9 

+48:00:15 

19.1 

3.297 

15 Apr 2005 

C5 

HIRESr 

10800 

13/10 







14 Jul 2005 

C5 

HIRESr 

10800 








04 Jun 2006 

C5 

HIRESr 

7200 


J1608+0715 


16:08:43.9 

+07:15:09 

16.60 

2.88 

11 Apr 2007 

Cl 

HIRESb 

9000 

11/26 

J1712+5755 


17:12:27.74 

+57:55:06 

17.46 

3.01 

09 Sep 2004 

Cl 

HIRESb 

3600 

xh2 







02 May 2005 

C5 

HIRESb 

3900 








19 Aug 2006 

Cl 

HIRESb 

3900 








20 Aug 2006 

Cl 

HIRESb 

3900 


J1733+5400 


17:33:52.23 

+54:00:30 

17.35 

3.43 

02 May 2005 

C5 

HIRESb 

5400 

X/30 







22 Aug 2007 

Cl 

HIRESr 

5400 

35/35 

J2123-0050 


21:23:29.46 

-00:50:53 

16.43 

2.26 

20 Aug 2006 

E3 

HIRESb 

21600 

30/67 

Q2126-1538 


21:29:12.2 

-15:38:41 

17.3 

3.27 

08 Sep 2004 

Cl 

HIRESb 

7200 

10/16 

LB2203-1833 


22:06:39.6 

-18:18:46 

18.4 

2.73 

09 Sep 2004 

Cl 

HIRESb 

5400 


Q2231-00 

LBQS 2231-0015 

22:34:08.8 

+00:00:02 

17.4 

3.025 

01 Nov 1995 

C5 

HIRESO 

14400 

30 

SDSS2303-0939 


23:03:01.5 

-09:39:31 

17.68 

3.455 

08 Nov 2005 

C5 

HIRESr 

7200 

25/29 

SDSS2315+1456 


23:15:43.6 

+14:56:06 

18.52 

3.377 

04 Jun 2006 

C5 

HIRESr 

4800 

16/11 







08 Nov 2005 

C5 

HIRESr 

4400 


SDSSJ2334-0908 


23:34:46.4 

-09:08:12 

18.03 

3.317 

18 Sep 2007 

Cl 

HIRESr 

14400 

28/X 

Q2355+0108 


23:58:08.6 

+01:25:06 

17.5 

3.40 

28 Oct 2005 

C5 

HIRESb 

7200 

21/30 







04 Jan 2006 

C5 

HIRESb 

6300 



^ Magnitude from the SDSS database (r-band) or as listed in the SIMBAD Astronomical Database (V-band). 

^ Decker employed. 

^ Median signal-to-noise per 3.0km s~^ pixel in the quasar continuum at ~ 5000A for the old HIRES detector (HIRESO), ~ 3400/4000A for 
HIRESb, and ~ 6 OOO/ 8 OOOA for HIRESr. An “X” indicates no wavelength coverage or that the S/N was compromised by an LLS. 
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TABLE 2 

JOURNAL OF MIKE OBSERVATIONS 


QSO 

Alt. Name 

RA 

(J2000) 

DEC 

(J2000) 

r/U^ 

(mag) 

^em 

Date 

(UT) 

Slit'’ 

(") 

Exp 

(s) 

S/Ng,„e 

(pix “) 

S/N,^ed 

(pix 1) 

QOOOl-2340 


00:03:45.0 

-23:23:46 

16.7 

2.262 

10 Sep 2005 

1.0 

3000 

3/27 

16/19 

J0103-3009 

LBQSOlOl-3025 

01:03:55.3 

-30:09:46 

17.6 

3.15 

02 Sep 2004 

1.0 

2400 

7/9 

8/10 







04 Sep 2004 

1.0 

2400 

4/7 

7/8 

SDSSJ0106+0048 


01:06:19.2 

+00:48:23.3 

19.03 

4.449 

26 Aug 2003 

1.0 

8000 

X/X 

8/9 

SDSSJ0124+0044 


01:24:03.8 

+00:44:32.7 

17.9 

3.834 

28 Aug 2003 

1.0 

8000 

X/4 

24/17 

SDSSJ0209-0005 


02:09:50.7 

-00:05:06 

16.9 

2.856 

10 Sep 2005 

1.0 

5700 

X/10 

11/9 

SDSSJ0244-0816 


02:44:47.8 

-08:16:06 

18.2 

4.068 

26 Aug 2003 

1.0 

5500 

X/2 

24/12 

HE0340-2612 


03:42:27.8 

-26:02:43 

17.4 

3.14 

02 Sep 2004 

1.0 

2400 

6/11 

12/19 







04 Sep 2004 

1.0 

2400 



SDSSJ0344-0653 


03:44:02.8 

-06:53:00 

18.64 

3.957 

28 Aug 2003 

1.0 

3000 

X/X 

12/8 

SDSS0912+0547 


09:12:10.35 

+05:47:42 

18.05 

3.248 

10 May 2004 

0.7 

3600 

2/3 

6/X 

SDSSJ0942+0422 


09:42:02.0 

+04:22:44 

17.18 

3.28 

03 Apr 2003 

0.7 

6000 

X/9 

14/XX 

HE0940-1050 


09:42:53.2 

-11:04:22 

16.6 

3.08 

08 May 2004 

1.0 

7200 

X/40 

35/X 

SDSSJ0949+0335 


09:49:32.3 

+03:35:31 

18.1 

4.05 

04 Apr 2003 

0.7 

4000 

4/4 

12/8 







05 Apr 2003 

0.7 

4000 

3/5 

11/8 







06 Apr 2003 

0.7 

4000 

2/5 

11/10 

SDSSJ1025+0452 


10:25:09.6 

+04:52:46 

18.0 

3.24 

05 Apr 2003 

0.7 

4000 

1/7 

12/9 







06 Apr 2003 

0.7 

4000 

1/7 

12/9 

SDSSJ1028-0046 


10:28:32.1 

-00:46:07 

17.94 

2.86 

12 May 2004 

1.0 

6484 

X/8 

15/6 

SDSSJ1032+0541 


10:32:49.9 

+05:41:18.3 

17.2 

2.829 

10 May 2004 

1.0 

7200 

4/13 

21/X 

SDSSJ1034+0358 


10:34:56.3 

+03:58:59 

17.9 

3.37 

04 Apr 2003 

0.7 

12000 

1/4 

15/12 

QllOO-264 


11:03:25.6 

-26:45:06 

16.02 

2.145 

16 May 2005 

1.0 

2000 

4/l7 

16/15 







18 May 2005 

1.0 

2000 

14/37 

12/9 

HS1104+0452 


11:07:08.4 

+04:36:18 

17.48 

2.66 

19 May 2005 

1.0 

2000 

2/9 

12/14 

SDSSJlllO+0244 


11:10:08.6 

+02:44:58 

18.3 

4.12 

05 Apr 2003 

0.7 

8000 

1/3 

10/7 

SDSSJ1136+0050 


11:36:21.0 

+00:50:21 

18.1 

3.43 

06 Apr 2003 

0.7 

8000 

2/7 

14/13 

SDSSJ1155+0530 


11:55:38.6 

+05:30:50 

18.1 

3.47 

10 May 2004 

1.0 

3600 

3/9 

12/10 

SDSSJ1201+0116 


12:01:44.4 

+01:16:11 

17.5 

3.23 

03 Apr 2003 

0.7 

8000 

1/5 

15/10 

LB1213+0922 


12:15:39.6 

+09:06:08 

18.26 

2.723 

13 May 2004 

0.7 

7200 

3/10 

10/10 

SDSSJ1249-0159 


12:49:57.2 

-01:59:28 

17.8 

3.64 

06 Apr 2003 

0.7 

8000 

1/13 

16/15 

SDSSJ1307+0422 


13:07:56.7 

+04:22:15 

18.0 

3.02 

09 May 2004 

0.7 

7200 

2/6 

10/10 

SDSSJ1337+0128 


13:37:57.9 

+02:18:20 

18.13 

3.33 

12 May 2004 

1.0 

6800 

X/10 

10/9 

SDSSJ1339+0548 


13:39:42.0 

+05:48:22 

17.8 

2.98 

10 May 2004 

1.0 

7200 

6^3 

14/12 

SDSSJ1402+0146 


14:02:48.1 

+01:46:34 

18.8 

4.16 

05 Apr 2003 

0.7 

8000 

1/2 

12/9 

SDSSJ1429-0145 

Q1426-0131 

14:29:03.0 

-01:45:18 

17.8 

3.42 

06 Apr 2003 

0.7 

8000 

2/12 

13/11 







17 May 2005 

1.0 

8000 



Q1456-1938 


14:56:50.0 

-19:38:53 

18.7 

3.16 

18 May 2005 

0.7 

7200 

5/10 

14/22 

SDSSJ1503+0419 


15:03:28.9 

+04:19:49 

18.1 

3.66 

09 May 2004 

0.7 

7200 

1/3 

8/7 

SDSSJ1558-0031 


15:58:10.2 

-00:31:20 

17.6 

2.83 

06 Apr 2003 

0.7 

6000 

1/6 

12/8 







10 May 2004 

1.0 

8000 

3/11 

18/16 

Q1559+0853 


16:02:22.6 

+08:45:36.5 

17.3 

2.269 

17 May 2005 

1.0 

4000 

4/17 

11/17 

SDSSJ1621-0042 


16:21:16.9 

-00:42:50 

17.4 

3.70 

03 Apr 2003 

0.7 

6000 

X/2 

11/10 







05 Apr 2003 

0.7 

3000 

X/7 

16/12 







06 Apr 2003 

0.7 

3600 

X/7 

18/14 







08 May 2004 

1.0 

3600 

3/11 


PKS2000-330 

Q2000-330 

20:03:24.1 

-32:51:44 

17.3 

3.77 

02 Sep 2004 

1.0 

4800 

12/35 

24/19 

B2050-359 


20:53:44.6 

-35:46:52 

17.7 

3.49 

18 May 2005 

1.0 

4800 

X/8 

10/10 

Q2126-1538 


21:29:12.2 

-15:38:41 

17.3 

3.27 

05 Sep 2004 

1.0 

4800 

9/25 

19/23 

HE2215-6206 


22:18:51.3 

-61:50:54 

17.5 

3.32 

02 Sep 2004 

1.0 

2400 

8/20 

16/14 







04 Sep 2004 

1.0 

4000 

7/17 

17/19 

SDSSJ2303-0939 


23:03:01.4 

-09:39:30 

17.68 

3.455 

28 Aug 2003 

1.0 

8000 

X/14 

23/20 

HE2314-3405 


23:16:43.2 

-33:49:12 

16.9 

2.96 

02 Sep 2004 

1.0 

2400 

2^1 

13/11 

SDSSJ2346-0016 


23:46:25.7 

-00:16:00 

17.77 

3.49 

27 Aug 2003 

1.0 

8000 

X/14 

21/26 







28 Aug 2003 

1.0 

3000 



HE2348-1444 


23:48:55.4 

-14:44:37 

16.7 

2.93 

02 Sep 2004 

1.0 

2400 

14/22 

30/33 

HE2355-5457 


23:58:33.4 

-54:40:42 

17.1 

2.94 

02 Sep 2004 

1.0 

2400 

17/7 

13/15 


^ Magnitude from the SDSS database (r-band) or as listed in the SIMBAD Astronomical Database (V-band). 

^ Slit width employed. For the blue (red) side, a l" slit yields a FWHM resolution of 10.7 (13.6) km s~^ for a source that fills the slit. 
^ Median signal-to-noise per 3.0km s“^ pixel in the quasar continuum at ~ 3400/4000A. An X designates no flux. 

^ Median signal-to-noise per 4.2km s~^ pixel in the quasar continuum at ~ 6 OOO/ 8 OOOA. An X designates no flux. 
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Relative Velocity (km/s) 


Relative Velocity (km/s) 


Relative Velocity (km/s) 


Fig. 1. — H I Lya (top row) and Ly/3 (bottom row) profiles for three LLSs with Nm > 10^® cm“^. The damping wings of Lyo are well- 
resolved in these SLLSs and the blue curves indicate for the best-estimates and uncertainty of the A^hi values. These are well constrained, 
even in poor S/N data. 


but not for a comprehensive LLS survey. 

2.3. Data Reduction 

The HIRES spectra were reduced with the HIRedujj^ 
software package, primarily as part of the KODIAQ 
project (Lehner et al. 2014). Briefly, each spectral im¬ 
age was Dias-subtracted, ftat-fielded with pixel flats, and 
wavelength-calibrated with corresponding ThAr frames. 
The echelle orders were traced using a traditional flat- 
field spectral image. The sky ba ckground was subtracte d 
with a b-spline algorithm (e.g. Bochanski et al. 2009), 
and the quasar flux was further traced and optimally ex- 
tracted with standard techniques. These spectra were 
flux normalized with a high-order Legendre polynomial 
and co-added after weighting by the median S/N of each 
order. This yields an individual, wavelength-calibrated 
spectrum for each night of observation in the vacuum 
and heliocentric frame. When possible, we then com¬ 
bined spectra from quasars observed on multiple nights 
with the same instrument configuration. 

Processing of the MIKE spectra used the MIRedux 
package now bundled within the XIDL software pack- 
asCl This pipeline uses algorithms similar to HIRedux. 
The primary difference is that the flux is estimated to¬ 
gether with the sky using a set of b-spline models which is 
demanded by the short 5" slits employed with MIKE. In 
addition, these data were fluxed prior to coaddition using 
the reduced spectrum of a spectrophotometric standard 


® http://www.ucolick.org/~xavier/HIRedux/index.html 
^ http://www.ucolick.org/~xavier/IDL/index.html 


(taken from the same night in most cases). Therefore, 
we provide both fluxed and normalized spectra from this 
instrument. 

Details on the data reduction of ESI a nd MagE spectra 


are provided in previous pub lications (jProchaska et al. 
2003b |Eumagalli et al.|[2013 ). 


All of the reduced and calibrated spectra are available 
on the project’s websit^ The Keck/HIRES spectra will 
also be provided in the nrst data release of the KODIAQ 
project (Lehner et al. 2015). 


3. LLS DEFINITION 

Before proceeding to analysis of the sample, we strictly 
define the Lyman limit system. There are three aspects 
to the definition: 


1. The velocity interval analyzed, which also corre¬ 
sponds to a finite redshift window. 

2. The Vhi value of the system. 

3. The spatial proximity of the LLS to other astro- 
physical objects (e.g. the background quasar or a 
foreground DLA). 

Of these three, the first has received the least attention 
by the community yet may be the most important. Es¬ 
tablishing a precise definition, however, is largely arbi¬ 
trary despite the fact that it may significantly impact 
the studies that follow. This includes the assessment 


http: / / www.ucolick.org/~xavier/HD-LLS/DRl 
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LLSJ094253.4-110425 z =2.917, logA^Hi = 17.50 [MIKE] 
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LLSJ094253.4-110425 z =2.917, logN^^ = ll.50 [MagE] 
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Fig. 2.— Plots of the Lyman limit (left) and the Lya (right) profile for the LLS at 2 : = 2.917 toward J0924-1104. We observed this 
system with both the MIKE (top) and MagE (bottom) spectrometers at LCO. The spectra show a break in the flux at A 5 ^ 3585A but 
residual flux down to A ps 3500A where the Lyman limit from a lower redshift absorber occurs. By modeling the flux decrement of the 
2 = 2.917 LLS, we establish a precise estimate of its total TVhi value. It is also evident that the H I Lya profile alone offers very little 
constraint on A^hi- 


of g as kinematics ([Prochaska fc: Wolfe 1997 1, me tallic- 
ity (Prochter et al. ' 2U1(J "Pumagalli et al. 2011a), and 
even the contribution ot LLS s to the cosniic mean free 
path ( Prochaska et al.||2014b ). In this paper (and future 
publications), we adopt an observationally-motivated ve¬ 
locity interval of ±500km s“^ centered on the peak op¬ 
tical depth of the H I Lyman series. Frequently, 

we estimate ^lls' from the peak optical depth of a low- 
ion, metal-line transition. An LLS, then, is all of the 
optically thick gas at |(5t!| < 500km s“^ from In 

practice, we have not simply summed the H I column 
densities of all Lya absorbers within this interval. In¬ 
stead, we have adopted the integrated Vhi estimate from 
the Lyman limit decrement or adopt Vhi from the anal¬ 
ysis of damping in the Lya profiles (see the next section 
for more detail). As an example of the latter, we treat 
the two absor bers at z = 3.1878 an d z = 3.1917 towards 
PKS2000-330 ( Prochter et al.|20l0 ) as a single LLS. Sim¬ 
ilarly, we sum metal-line absorption identified within the 
interval although it rarely is detected in intervals that 
exceeds 200km s“^. Moreover, this window was adjusted 
further to exclude absorption from unrelated (e.g higher 
or lower redshift) systems. While this is an observation- 
ally driven definition, we note that it should also capture 
even the largest peculiar motions within dark matter ha¬ 
los at z ^ 2. 

With |(5z;| < 500km s“^ as the first criterion, we de¬ 
fine an LLS as any combination of systems with Vhi > 


lO^"^'^ cm“^ within that interval; this yields an integrated 
optical depth at the Lyman limit Tgi 2 >1. In prac¬ 
tice, we distinguish the LLS from DLAs by requiring 
that Vhi < 10^°-^cm“^. Systems with 10^® < Am < 
10^^'^ cm“^ are referred to as partial LLSs or pLLSs, and 
are excluded from analysis in this manuscript. Last, we 
refer to an LLS within 3000 km s“^ o f the background 


quasar as a proximate LLS or PLLS (Prochaska et al. 


2008b I. There are 5 PLLSs within our sample satisfying 


this definition, all with velocity separations of at least 
2000km s“^ from the reported quasar redshifts. Alto¬ 
gether, we present measurements for 157 LLSs at red- 
shifts zlls = 1-76 — 4.39 and with Whi = 10^^'^ — 
1020.25 cm“2 Here and in future papers we refer to 
this dataset as the high-dispersion LLS sample (HD-LLS 
Sample). We will augment this sample in the years to 
follow via our web site. 


4. Afni ANALYSIS 

Although the continuum opacity of the Lyman limit 
generates an unambiguous signature in a quasar spec¬ 
trum, it is generally challenging to precisely estimate the 
H I column density Vhi for a given LLS system. This 
follows simply from the fact that exp(—Tgi 2 ) <C 1 for 
Vhi > 10^® cm“^ and all of the H I Lyman series lines 
are on the saturated portion of the curve-of-growth for 
Vhi < 10^®cm“^. Furthermore, the damping of Lya 
is difficult to measure for Nm <C 10^°cm“^, especially 
in low S/N spectra or at z > 3 where IGM blending is 
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LLS J235057.87-005209.9 z =2.929, logA^ni =(17-40,18.50,18.90) 
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Fig. 3.— Two examples of LLSs with a strong Lyman limit break (r 9 i 2 > 1) yet weak or absent damping of H I Lya. In these examples, 
we only estimate bounds on the A^hi values which can span an order-of-magnitude uncertainty. 


substantial. 

Our approach to identifying LLS and estimating their 
A^hi values involved two relatively distinct procedures. 
For LLS with large A^hi values (> 10^®cm“^), we 
searched each spectrum for absorption features with large 
equivalent widths Wx characteristic of a damped Lya line 
(i.e. Wx lA). We then considered whether these can¬ 
didates could be related to a broad absorption line (BAL) 
system associated to the background quasar or Ly/? asso¬ 
ciated to a higher redshift DLA. Any such coincidences 
were eliminated. For the remaining candidates, we per¬ 
formed an analysis of the Lya prohle by overplotting a 
series of Voigt profiles with Vhi > 10^®'^ cm“^, adjusting 
the local continuum by-eye as warranted. When low-ion 
metal absorption was detected near the approximate cen¬ 
troid of Lya, we centered the model to its peak optical 
depth and refined the Vri value accordingly. We did 
not, however, require the positive detection of metal-line 
absorption. In all cases, the Doppler parameter of the 
model Lya line was set to 30 km s“^. 

For cases were the S/N was deemed sufficient and line¬ 
blending not too severe, we estimated (by visual inspec¬ 
tion) a ‘best’ A^hi value and corresponding Icr uncertain¬ 
ties. Although this procedure is rife with human inter¬ 
action, we maintain that it offers the most robust as¬ 
sessment (to date) for A^hi estimation. This is because 
the dominant uncertainties are systematic (e.g. contin¬ 
uum placement and line-blending), which are difficult to 
estimate statistically. Figure shows three examples of 
LLSs with damped Lya lines giving precisely estimated 
A^hi values. Such systems are commonly referred to as 


super LLSs (SLLSs) or sub-DLAs. 

We provide the adopted A^hi values and error estimates 
of our SLLS sample (99 systems with A^hi A 10^® cm“^) 
in Table These represent roughly 2/3 of the total 
HD-LLS Sample. This high fraction occurs because we 
only require coverage of H I Lya to identify and ana¬ 
lyze these SLLS. This implies a much larger survey path- 
length than for the lower A^hi LLS. In addition, one may 
identify and analyze multiple SLLSs along a given sight- 
line whereas one is restricted to a single LLS when the 
Lyman Limit is central to the analysis. Because these 
SLLSs tend to span nearly the entire ±500km s“^ win¬ 
dow that defines an LLS, it is possible that there is ad¬ 
ditional, optically thick gas not included in our A^hi esti¬ 
mate. This will be rare, however, and the underestimate 
of A^hi should generally be much less than 10%. 

In parallel with the search for LLSs having strong Lya 
lines, we inspected each spectrum for a Lyman limit 
break. For those the LLSs that exhibit non-negligible 
flux at the Lyman limit, i.e. tqi 2 < 3, a precise A^hi 
estimation may be recovered. In practice, such analy¬ 
sis is hampered by poor sky subtraction and associated 
IGM absorption that stochastically reduces the quasar 
flux through the spectral region near the Lyman limit 
and affects continuum placement. In the following, we 
have been conservative regarding the systems with A^hi 
measurements from the Lyman limit flux decrement. We 
are currently acquiring additional, low-dispersion spectra 
to confirm the flux at the Lyman limit for a set of the 
HD-LLS Sample. Figure shows one example of a pLLS 
observed with both the MIKE and MagE spectrometers. 


































































































































































































































































































TABLE 3 

a^hi estimates for the hd-lls sample 


Quasar 

peak 

AfloW 

;u-high 

»7-adopt 

/^5hi 

J1608+0715 

1.7626 

19.10 

19.70 

l9.40/g;30 

1 

J0953+5230 

1.7678 

20.00 

20.20 

20 .l0/g;(° 

1 

J0927+5621 

1.7749 

18.90 

19.10 

l9.00/g;(o 

1 

J1509+1113 

1.8210 

18.00 

19.00 

l8.50/g;5° 

2 

J101939.15-I-524627 

1.8339 

18.80 

19.40 

I9.l0/g;®° 

1 

QllOO-264 

1.8389 

19.25 

19.55 

l9.40/g;(5 

1 

J1159-0032 

1.9044 

19.90 

20.20 

20.05/®;)® 

1 

Q0201-I-36 

1.9548 

19.90 

20.30 

20 .10/“;®S 

1 

J0828+0858 

2.0438 

19.80 

20.00 

19.90/g;(° 

1 

J2123-0050 

2.0593 

19.10 

19.40 

19.25/®;)® 

1 

Q1456-1938 

2.1701 

19.55 

19.95 

19 75+0-20 
ly. < o_Q 20 

1 

J034024.57-051909 

2.1736 

19.15 

19.55 

19.35/g;®° 

1 

QOOOl-2340 

2.1871 

19.50 

19.80 

19.65/g;(® 

1 

SDSS1307-I-0422 

2.2499 

19.85 

20.15 

20 .00/g;(® 

1 

J1712+5755 

2.3148 

20.05 

20.35 

20.20/g;(5 

1 

Q2053-3546 

2.3320 

18.75 

19.25 

l9.00/g;®5 

1 

Q2053-3546 

2.3502 

19.35 

19.85 

19.60/°;25 

1 

Q1456-1938 

2.3512 

19.35 

19.75 

19 55+0-20 
iy.oo_Q 20 

1 

J1131+6044 

2.3620 

19.90 

20.20 

20.05/g;(5 

1 

Q1206-I-1155 

2.3630 

20.05 

20.45 

20.25/“;®S 

1 

HE2314-3405 

2.3860 

18.80 

19.20 

l9.00/g;20 

1 

Q0301-005 

2.4432 

19.75 

20.05 

19.90/g;(5 

1 

HS1345+2832 

2.4477 

19.70 

20.00 

19.85/g;(5 

1 

J1035+5440 

2.4570 

19.40 

19.90 

19.65/g;®5 

1 

Q1337-I-11 

2.5080 

20.00 

20.30 

20.15/g;(® 

1 

SDSS0912-I-0547 

2.5220 

19.15 

19.55 

19.35/g;®° 

1 

SDSS0209-0005 

2.5228 

18.90 

19.20 

19-0+2+ 

1 

LB1213+0922 

2.5230 

20.00 

20.40 

20 .20/“;®S 

1 

Q0207-003 

2.5231 

18.80 

19.20 

19.00/“;2o 

1 

Q0207-003 

2.5466 

17.60 

18.60 

18-10/2:22 

2 

HS1104+0452 

2.6014 

19.70 

20.10 

19.90/2:22 

1 

J2234+0057 

2.6040 

19.25 

19.75 

19.50/2:22 

1 

J115659.59-I-551308.1 

2.6159 

18.80 

19.30 

19 10+0-30 
iy.iu_o 30 

1 

SDSS1558-0031 

2.6300 

19.40 

19.75 

iy.DU_Q 20 

1 

SDSS0157-0106 

2.6313 

19.25 

19.65 

19.45/2:22 

1 

Q2126-158 

2.6380 

19.10 

19.40 

19.25/2:22 

1 

Q1455-I-123 

2.6481 

17.30 

19.40 

18-35/1:22 

2 

LBQS2231-0015 

2.6520 

18.70 

19.30 

iy.iu_o 40 

1 

SDSS0121-I-1448 

2.6623 

19.05 

19.40 

19.25/2:22 

1 

SDSSJ0915-I-0549 

2.6631 

17.50 

18.90 

18.20/2:12 

2 

SDSSJ2319-1040 

2.6750 

19.30 

19.60 

19.45/2:12 

1 

Q0201-I-36 

2.6900 

17.50 

18.80 

18.50/2:22 

2 

LB2203-1833 

2.6981 

19.85 

20.15 

20 .00/2:12 

1 

SDSSJ1551-I-0908 

2.7000 

17.30 

17.70 

17.50/2:22 

1 

HS1200+1539 

2.7080 

17.60 

18.90 

18-30/2:12 

2 

Q1508-I-087 

2.7219 

19.00 

19.80 

19.40/2:12 

1 

PMNJ1837-5848 

2.7289 

17.50 

18.70 

18-10/2:22 

2 


Note. —All column densities are logio. The flag in the final column 
indicates the quality of the measurement. A flgui = 1 corresponds 
to a more precisely measured value and one may assume a Gaussian 
PDF with the errors reported taken as la uncertainties. A flgui = 2 
corresponds to a less precisely measured value, and we recommend one 
adopt a uniform prior for Nm within the error interval reported. See 
text for further details. 

Redshift estimates for the peak H I opacity from metals and Lyman 
series absorption. 



Fig. 4.— (left) Histogram of A^hi values for the HD-LLS Sample. 
The ligher bars indicate Afjji values with large uncertainty (Icr > 
0.4dex). Nearly two-thirds of the sample have Nhi > 10^®cm“^ 
which is a consequence of the spectral coverage required to estimate 
Nm for an LLS (see text), (right) Scatter plot of A^hi vs. zlls fo'" 
the sample. Again, the lighter points indicate LLSs with poorer 
constraints on the AIhi values (flg=2 in Table |^. Systems with 
Nm < cm“® all have zlls > 2.6 as coverage of the Lyman 
limit is required for the H I analysis. 


The flux decrement is obvious and one also appreciates 
the value of higher spectral resolution (with high S/N) 
for resolving IGM absorption. 

For the remainder of the systems identified on the 
basis of a Lyman limit break, we adopt conservative 
bounds (i.e. upper and lower limits) to the iVni values. 
These are based primarily on analysis of the Lyct line 
and the flux at the Lyman limit. The absence of strong 
damping in the former provides a strict upper limit to 
A^hi while the latter sets a firm lower limit. These 
bounds are provided in Table and Figure shows 
two examples of these ‘ambiguous’ cases. In practice, 
the bounds are often an order-of-magnitude apart, e.g. 
lO^^'^cm”^ < iVni < 10^® ®cm“^. Furthermore, it is 
difficult to estimate the probability distribution function 
(PDF) of A^hi within these bounds. One should not, for 
example, assume a Gaussian PDF centered within the 
bounds with a dispersion of half the interval. In fact, 
we expect that the PDF is much closer to uniform, i.e. 
equal probability for any A^hi value within the bounds. 
This expectation is motivated by current estimations of 
the A^hi frequency distribution /(A^hi,A') which argue 
for a uniform distribution of A^hi values for randomly se- 
lected systems with A^ht ~ 10^®cm“^ (Prochaska et al.l 


2010 

O’Meara et al. 

2013 

adopting a unilorm J 



As a cross-check on the analysis, 50 of the sightlines 
were re-analyzed by a second author to identify LLSs and 
estimate their Ahi values. With two exceptions, the val¬ 
ues between the two authors agree within the estimated 
uncertainty and we identify no obvious systematic bia^ 
These two exceptions have > 0.5 dex difference due to dn- 
fering definitions used by the two authors and we have 
adopted the values corresponding to the strict definition 
provided in § |^ This exercise confirms that the uncer¬ 
tainties are dominated by systematic effects, not S/N nor 
the analysis procedures. 

Table lists the adopted A^hi value, errors on this 
value, the bounds on A^hi, and a flag indicating whether 
one would assume a normal or uniform PDF. Figure 


® In fact the formal reduced x® for the comparison is significantly 
less than unity, but this is because the estimated uncertainties in¬ 
clude systematic error and because each author analyzed the same 
data. 
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Relative Velocity (km/s) 


Relative Velocity (km/s) 


Relative Velocity (km/s) 


Fig. 5. — Si II and Si IV transitions for three SLLS representative of the full HD-LLS Sample. Note the large diversity in metal-line 
strength despite the comparable Nm values. These systems have a tendency to sho-w both low and high-ion absorption indicative of 
partially ionized gas. The gray dotted line in each panel indicates the estimated Icr error array. 


TABLE 4 

IONIC COLUMN SUMMARY FOR Si AND C 


Quasar 

^ahs 

A^hi 

N{C+) 

a{N{C+)) 

7V(C°+) 

a{N{C^+)) 

7V(Si+) 

^(lV(Si+)) 

Ar(Si°+) 

<T(7V(Si3+)) 

J1608+0715 

1.7626 

19.40l°j° 





15.80 

-9.99 



J0953+5230 

1.7678 

20 .10t°;;° 

15.44 

+9.99 

15.21 

+9.99 

15.67 

0.01 

14.57 

+9.99 

J0927+5621 

1.7749 

19.001°;;° 

15.40 

+9.99 

15.40 

+9.99 

15.58 

0.02 

14.84 

+9.99 

J1509+1113 

1.8210 

18.501°;5° 



14.83 

+9.99 

14.21 

0.04 

14.17 

+9.99 

J101939.15+524627 

1.8339 

19.101°;°° 



14.93 

+9.99 

15.32 

0.03 

14.14 

+9.99 

QllOO-264 

1.8389 

19.401°;;° 



14.24 

0.00 

13.96 

0.01 

13.83 

0.00 

J1159-0032 

1.9044 

20.051°;;° 

15.38 

+9.99 

15.22 

+9.99 

15.14 

0.10 

14.54 

+9.99 

Q0201+36 

1.9548 

20 .101°;"° 





15.11 

0.09 



J0828+0858 

2.0438 

19.901°;;° 

15.14 

+9.99 

14.89 

+9.99 

15.25 

0.10 

14.44 

+9.99 

J2123-0050 

2.0593 

19.251°;;° 

15.11 

+9.99 

14.60 

+9.99 

14.60 

0.04 

13.96 

0.00 

Q1456-1938 

2.1701 

IQ ■7C+0.20 
iy. t o_Q 20 





14.84 

-9.99 



J034024.57-051909 

2.1736 

19.351°;^° 

14.40 

+9.99 

13.86 

0.02 

13.84 

0.02 

13.39 

0.02 

QOOOl-2340 

2.1871 

19.651°;;° 

14.45 

+9.99 

14.26 

0.01 

13.75 

0.03 

13.74 

0.01 

SDSS1307+0422 

2.2499 

20 .001°;;° 



14.22 

0.03 

14.25 

+9.99 



J1712+5755 

2.3148 

20 .201°;;° 



13.36 

0.04 

14.08 

0.01 




Note. — [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 
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Relative Velocity (km/s) Relative Velocity (km/s) 


Relative Velocity (km/s) 


Fig. 6. — C II and C IV transitions for LLSs with low A^hi values. Unlike the SLLSs from Figure!^ 
that is dominated by high-ions. There only a few cases of LLSs with A^hi < 10^® cm~^ and a positiv^ 


these LLSs have metal absorption 
ow-ion detection. 


TABLE 5 

IONIC COLUMN SUMMARY FOR Al, Fe, and O 


Quasar 

•^abs 

A^hi 

N{0°) 

<7(7V(o0)) 

Af(Al+) 

a{N{Al+)) 

Af(Al++) 

cr{N{Al++)) 

Ar(Fe+) 

a{N(Fe+)) 

J1608+0715 

1.7626 

19.40t°;|° 





13.53 

0.00 



J0953+5230 

1.7678 

20.101°;;° 

15.68 

+9.99 

13.96 

+9.99 

13.86 

0.01 

14.99 

0.10 

J0927+5621 

1.7749 

19.001°;;° 

15.63 

+9.99 

13.92 

+9.99 

14.05 

0.01 

15.28 

0.13 

J1509+1113 

1.8210 

18.50l°'j° 



13.12 

+9.99 

13.04 

0.05 

13.76 

0.11 

J101939.15+524627 

1.8339 

19.101°;°° 



13.34 

+9.99 

13.62 

0.02 

14.19 

0.02 

QllOO-264 

1.8389 

19.401°;;° 



12.79 

0.01 

12.31 

0.10 

13.42 

0.01 

J1159-0032 

1.9044 

20.051°;;° 

15.66 

+9.99 

13.98 

+9.99 

13.82 

0.01 



Q0201+36 

1.9548 

20.101°;°° 



13.77 

+9.99 

13.61 

0.01 



J0828+0858 

2.0438 

19.901°;;° 

15.49 

+9.99 



13.59 

0.02 

14.89 

0.04 

J2123-0050 

2.0593 

19.251°;;° 



13.44 

+9.99 

13.15 

0.01 

14.39 

0.00 

Q1456-1938 

2.1701 

IQ ■7C+0.20 
iy. 1 o_Q 20 



13.38 

+9.99 

12.99 

0.04 

14.26 

0.01 

J034024.57-051909 

2.1736 

19.351°;°° 

14.56 

+9.99 

12.65 

0.03 

12.49 

0.14 



QOOOl-2340 

2.1871 

19.651°;;° 

14.16 

0.04 

13.00 

+9.99 

12.40 

-9.99 

13.11 

0.03 

SDSS1307+0422 

2.2499 

20.001°;;° 



13.04 

+9.99 

12.80 

0.07 

14.18 

0.04 

J1712+5755 

2.3148 

20.201°;;° 



12.56 

0.02 

12.39 

0.06 

13.65 

0.03 

Note. — [The complete version 

of this table 

is in the electronic edition of the Journal. The printed edition contains only 

a sample.] 
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shows a histogram of the adopted A^hi values and a scat¬ 
ter plot against zlls- It is evident that the HD-LLS 
Sample is weighted towards higher A^hi values and z ^ 3. 

5. IONIC COLUMN DENSITIES 

For each of the HD-LLS Sample, we inspected the 
spectra for associated metal-line absorption. Empha¬ 
sis was placed on transitions with observed wavelengths 
redward of the Lya forest. A velocity interval was 
estimated for the column density measurements based 
on the cohort of transitions detected. Velocity plots 
were generated and inspected to search for line-blending. 
Severely blended lines were eliminated from analysis and 
intermediate/weak cases were measured but reported 
as upper limits. All of these assignments were vet¬ 
ted by JXP, MF, and JMO. Figures and show the 
Si II 1526/Si IV 1393 transitions for three representa¬ 
tive SLLSs and the C II 1334/C IV 1548 transitions for 
three LLSs with Vhi ^ 10^^'^ cm“^. These data indicate 
a great diversity of line-strengths for these transitions 
within the SLLS sample. We also conclude that metal- 
absorption is dominated by high-ions in the lower Vhi 
systems. 

Column densities were mea sured using the apparent 
optical depth method (AODM; Savage &: Sembach|1991| ) 
which gives accurate results for unsaturated line profiles. 
On the latter point, the echelle data (MIKE,HIRES) 
have sufficiently high resolution to directly assess line- 
saturation, i.e. only profiles with minimum normalized 
flux /min less than 0.1 may be saturated. For the echel- 
lette data (MagE, ESI), however, l ine-saturation is a con¬ 
cern (e.g. Prochaska et al. 2003bI. In general, we have 
proceeded conservatively by treating most lines as satu¬ 
rated when /min < 0.5. For many of the ions analyzed 
in these LLSs, we observe multiple transitions with dif¬ 
fering oscillator strengths and have further assessed line- 
saturation from the cohort of measurements. 

Uncertainties were estimated from standard propaga¬ 
tion of error, which does not include error from contin¬ 
uum placement. To be conservative, we adopt a mini¬ 
mum uncertainty of 0.05 dex to the measurements from 
a given transition. When multiple transitions from the 
same ion were measured (e.g. Si II 1304 and Si II 1526 
for Si^), we calculate the ionic column density from the 
weighted mean. Otherwise, we adopt the measurement 
from the single transition or a limit from the cohort em¬ 
phasizing positive detections. 

A complete set of tables and figures for the metal-line 
transitions analyzed for each LLS are given online. Ta¬ 
bles 0 and summarize the results for A1+, A1++, Fe^, 
C’*', CU, Si+, and Si’*'^. A listing of all the measure¬ 
ments from this manuscript is provided in the Appendix. 
Figures and show the column density distributions 
for a set of Al, Fe, Si, C, and O atoms/ions as a func¬ 
tion of the LLS A/ji value. Not surprisingly, the lower 
ionization states show an obvious correlatiorrD with H I 
column density although there is a large scatter at all val¬ 
ues. The near absence of positive detections for O I (i.e. 
V(0°) < 10^“^ cm“^) at low Vpi is also notable. This 
suggests a rarity of high metallicity gas in systems with 
t '912 < 10. The high ions are also positively correlated 


Taking limits as values, all of these ions have a Spearman 
rank test probability of less than 0.0001. 



Fig. 7.— Scatter plot of Si and C ionic column densities for 
the HD-LLS Sample. Circles indicate measured values; their un¬ 
certainties are generally less than 0.1 dex. Triangles indicate lim¬ 
its to the values with the open symbols indicating upper limits. 
Lighter points mark LLSs with a poorly constrained ^hi value. 
Gray boxes encompass 50-percent of the measurements in three 
logarithmic Nm intervals: [17.3, 18.0), [18.0, 19.0), [19.0, 20.3). 
At all column densities, there is a large dispersion in the mea¬ 
surements. Nevertheless, the low-ions (C"^, Si"^) exhibit a strong, 
positive correlation with Nm value. A Spearman rank test rules 
out the null hypothesis at > 99.99% c.l. 



Fig. 8 .— Same as Figure [7| but for four additional ions. 

with the neutral column but with yet larger scatter and 
much smaller correlation coefficients. 

6 . RESULTS 

In the following, we present a set of results derived 
from the column density measurements of the previous 
sections. For this manuscript, we restrict the analysis 
to an empirical investigation. Future studies will intro¬ 
duce additional models and analysis (e.g. photoionization 
modeling) to interpret the data. We further restrict the 
discussion to ionic abundances and defer the analysis of 
kinematics to future work. 
















12 




Fig. 9. — Scatter plots of four ionic ratios that diagnose the 
ionization state of the LLSs. Gray boxes encompass 50-percent of 
the measurements in two logarithmic Afjji intervals: [17.3, 19.0), 
[19.0, 20.3). The C, Si, and A1 ratios show strong evidence that the 
gas is more highly ionized at low Nm values. Similarly, the set of 
Si+ZO'^ values exceeding —0.7 dex are indicative of a highly ionized 
gas. We further emphasize that the log(Si+/Si+®) ~ —0.5 dex 
values at high A^jjl suggests that this gas is also partially ionized. 

6.1. Ionization State 

As noted above, a full treatment of the ionization state 
of the gas including the comparison to models will be pre¬ 
sented in a future manuscript. We may, however, explore 
the state of the gas empirically through the examination 
of ionic ratios that are sensitive to the ionization state of 
the gas. Figure [^presents four such ionic ratios against 
A^hi- These primarily compare ions of the same element 
(e.g. to eliminate offsets due to differing intrin¬ 

sic chemical abundances (i.e. varying abundance ratios). 
In this analysis, we have taken the integrated column 
density across the entire LLS. While there is evidence 
for variations in these ratios with in individual compo- 
nents, these tend to be small (e.g. Prochter et al. 2010 
Figure Therefore, the trends apparent in 1'igure 
reflect the gross properties of the LLS sample. 

All of the C+VC+, Si+3/Si+, and A1++/A1+ ratios 
exhibit a strong anti-correlation with Ahi indicating an 
increasing neutral fraction with increasing H I opacity. 
Taking limits as values, the Spearman rank test yields 
a probability of less than 10“^ for the null hypothesis, 
in each case. For all of these ions, the upper ionization 
state is dominant for A^hi [5 10^®'^ cm“^ and vice-versa 
for higher A^hi values. We emphasize, however, that even 
at A^hi ~ 10^° cm“^ the observed ratios are frequently 
large, e.g. log(Si+®/Si''") « —0.5 dex. This suggests that 
the gas is predominantly ionized even at these larger total 
H 1 opacities. 

This inference is further supported by the Si+/0°. 
Ignoring differential depletion, which we expect to be 
modest in LLSs, the Si+/0° ratio should trend towards 
the solar abundance ratio (esi/eo = —1.2 dex) in a 
neutral gas given that Si and O are both produced in 
massive stars and are observed to trace each other in 
astrophysical systems (e.g. stellar atmospheres). We 
identify, however, a significant sample of systems with 


5 . 

o 


-1 


-2 


-3 



18 19 

log Vhi 


20 


18 19 

log Vhi 


Fig. 10.— Scatter plots of low-ion column densities relative 
to H I, normalized to the solar abundance and plotted 

against the LLS A^hi value. Gray boxes encompass 50-percent of 
the measurements in two logarithmic A^hi intervals: [17.3, 19.0), 
[19.0, 20.3). If ionization corrections are small, provides 

an estimate of the logarithmic metal abundance relative to so¬ 
lar. The measurements appear to indicate a declining trend of 
gas metallicity with increasing A/’hi- We argue, however, that this 
apparent trend is driven by ionization effects and the set of up¬ 
per/lower limits at low/high TVhi values. Furthermore, given the 
large scatter at all TVhi values, it will be challenging to establish 
any trend between enrichment and TVhi value in the LLSs. 

Nm « 10^®'®cm“^ that have log(Si+/0°) > —Idex. Be¬ 
cause the majority of ionization processes (e.g. photoion- 
izati on, collisional ionization) pre dict Si'''/0° > Si/0 
(e.g. Prochaska & Hennawi 20091, these measurements 
offer further evidence that LLSs are highly ionized. 

6.2. Metallicity 

A principal diagnostic of the LLSs is the gas metallic¬ 
ity, i.e. the enrichment of the gas in heavy elements. This 
quantity is generally characterized relative to the chemi¬ 
cal abundances observed for the Sun. For the f ollowing, 
we adopt th e solar abundance scale compiled by Asplund| 
et al. (2009), taking meteoritic values when possible. 

Because the LLSs are significantly ionized, the ob¬ 
served ionic abundances reflect only a fraction of the total 
abundances of Si, O, H, etc. Therefore, a full treatment 
requires ionization modeling. We may, however, offer in¬ 
sight into the problem by examining several ions relative 
to H*^. To minimize ionization corrections, one restricts 
the analysis to ionization states dominant in a highly 
optically thick (i.e. neutral) medium. 

The results for four low-ions are presented in Fig¬ 
ure!^ normalized to the solar abundance. We have in¬ 
troduced here a new quantity and notation: {XYY-^} = 
log(V(X*)/ZV(Yj)) — ex + ey, where ex is the solar abun¬ 
dance on the logarithmic scale for element X. This quan¬ 
tity explicitly ignores ionization corrections and should 
not be considered a proper estimate of the chemical 
abundance ratio, traditionally expressed as [X/Y]. In the 
cases where ionization corrections are negligible, how¬ 
ever, {XyH®} = [X/H] and this quantity represents the 
metallicity. 

A cursory inspection of the plots suggests a signifi- 
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Fig. 11.— Estimations of the oxygen metallicity in the LLSs 
where we have assumed that [0/H]={0'^/H*'}, i.e. that ionization 
corrections are small for this ionic ratio (see the text). Gray boxes 
encompass 50-percent of the measurements in three logarithmic 
A^hi intervals: [17.3, 18.0), [18.0 , 19 .0), [19.0, 20.3). Similar to the 
results for other low-ions (FigurellM, there is an apparent decline in 
[O/H) with Niii. This apparent trend, however, is primarily driven 
by the high incidence of upper/lower limits at low/high Nm values. 
In fact, the underlying trend may well be the opposite. Compare 
the set of low [O/H) values and upper limits at Afni ~ 10^® cm~^ 
with the large set of lower limits at Nm ^ 10^®'^cm“^. 

cant decline in metal content with increasing JVhi- This 
apparent anti-correlation, however, is driven by at least 
two factors. First, larger A^hi implies larger metal col¬ 
umn densities such that the transitions saturate yielding 
a preponderance of lower limits. By the same token, at 
low /Vhi values the transitions are often undetected yield¬ 
ing upper limits to the ionic ratios. Second, we have ar¬ 
gued from Figure that the gas is increasingly ionized 
with decreasing i\%i. For Si+, 0+, and A1+, the ion- 
ization corrections for |XVH°| are lik ely negative (e.g. 
Prochaska|1999 Fumagalli et al. 2011a), and would lower 
the metallicity one inters from such ratios. We believe 
these factors dominate the trends apparent in Figure [lOl 
In fact, it is even possible that the true distribution ex¬ 
hibits the opposite trend. Figure [TI] shows [O/H] against 
A^hi for the LLSs where we have assumed no ioniza¬ 
tion corrections, i.e. [O/H] = {0°/H°}. This approx¬ 
imation is justified by the fact that 0° and H° have 
very similar ionization potentials and their neutral states 
are coupled by charge-exchange reactions. This assump¬ 
tion may break down at low A ^ht values in the pres¬ 


ence of a hard radi ation field (Sofia & Jenkins] 1998 
Prochter et al. 2010), but the corrections are still likely 
to be modest (several tenths dex). Unfortunately, the 
measurements are dominated by limits: upper limits at 
A^hi < 10^®'® cm“^ and lower limits at A^hi > 10^® cm“^. 
Nevertheless, the data require that [O/H] > —1.7 for 
the SLLSs and indicate [O/H] < —1.3 dex for LLSs with 
A^hi ~ 10^®cm“^. We tentatively infer that the me¬ 
dian metallicity is approximately fiat with A^jji and pos¬ 
sibly increasing; more strictly, we rule out a steeply de¬ 
clining O/H metallicity with increasing A^hi- A similar 
conclusion may be drawn from the {Si~*’/H°} measure- 
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Fig. 12.— Comparison of {X*/H®} measurements for C"*" and Si“^ 
against absorption redshift. Left-hand panels are for the LLSs with 
A^hi < 10^^ cm“^ and the right-hand panels are for the SLLS pop¬ 
ulation. All of the data show evidence for a declining enrichment 
with increasing redshift, although this assertion is statistically sig¬ 
nificant (> 99% c.l.) only for the SLLS sample. The absence of 
low values at z ^ 2 implies a reduced incidence of near-pristine gas 
with high H I columns at that epoch. 


ments which scatter less from line-saturation. The LLSs 
with Nhi « 10^® cm“^ show very few positive detec¬ 
tions and have a median {Si^/H°} < — Idex. In con¬ 
trast, the LLSs at Nhi ~ 10^® ® cm“^ frequently exhibit 
{Si+/H0} > -Idex. 


Another result apparent from Figure 10 is the large 
dispersion in measurements at every Nhi value. This 
is most notable for Si’'" which has multiple transitions 
that permit measurements of the column density over 
a larger dynamic range. At the largest Nhi values, 
the values/limits of {Si’'’/H*^} span nearly four orders 
of magnitude! And although the measurements for 
LLSs with Nhi ~ ® — 10^® cm“^ include many up¬ 

per limits, one identifies values and upper limits with 
{Si’'’/H®} > —0.5 dex together with upper limits having 
{Si+/H0} < —2. Clearly, any underlying trend of en¬ 
richment with Nhi will be diluted by the large intrinsic 
scatter within the LLSs. One may even argue that if 
such a dispersion is indicative of multiple astrophysical 
systems, then defining a mean of the LLS population has 
limited scientific value. 

Despite the large dispersion, we emphasize that very 
few of the LLS in the HD-LLS Sample are “metal-free”, 
i.e. exhibiting no metal-line absorption and therefore con¬ 
sistent with primordial abundances. Of the tills LLSs, 
only 25 have no low-ion detections outside the Lya for¬ 
est and 18 of these exhibit a positive detection in a 
higher-ion. For the other 7, one has been previously been 


identifi ed as consistent with primordial (Fumagalli et al. 
2011a). The remainder have a diversity of S/iN and spec¬ 


tral coverage and therefore are generally less sensitive to 
measuring a low metallicity. Several will be examined 
in greater detail in a future manuscript. Nevertheless, 
we may conclude that the incidence of very low metal¬ 
licity gas (< 1/1000 solar) is rare in the LLS popula- 
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tion (< 5%). Furthermore, none of the 82LLSs with 
A^hi > 10^®'^ cm“^ are metal-freej^ By z ~ 3, gas that 
is dense enough to exhibit a very high Lyman limit opac¬ 
ity has previously been polluted by heavy elements. 

At the opposite end of the enrichment distribution, we 


identify 13 systems with a positive {Si''’/H°} measure¬ 
ment that exceeds 0 dex. This includes four extreme ex¬ 
amples with {Si'''/H°} > -1-0.5dex. Because these four 
LLSs also have Nm > 10^® cm“^ we expect that cor¬ 


rectio ns for ionization are modest (see Prochaska et al. 


20061 and that these are truly super-solar abundances. 


The others, however, have uncertainties consistent with 
the gas being sub-solar even before accounting for ioniza¬ 
tion. We conclude, subject to additional future analysis, 
that super- sola r enrichment is also rare in the LLSs. 

In Figure 12 we examine {Si'''/H°} and {C^/H^} val¬ 


ues as a function of redshift, splitting the LLS sample at 
A^hi = 10^® cm“^. The values for the lower A^hi systems 
suggest a declining trend with increasing redshift, e.g. in 
contrast to the lower redshift systems, none of the z > 3.5 
LLSs have a positive detection of {C+/H°} > —0.5 dex. 
Even if we restrict analysis to positive detections, how¬ 
ever, an anti-correlation is not statistically significant. 

Turning to the SLLS population, the {XYH°} distribu¬ 
tions show obvious trends with redshift (limits not with¬ 
standing) . Treating all of the positive detections at their 
plotted values, a Spearman’s rank correlation test rules 
out the null hypothesis at > 99% c.l. We interpret this 
anti-correlation as lower average enrichment within the 
SLLS at higher redshift. This conclusion relies on the as¬ 
sumption that ionization corrections will not evolve sig¬ 
nificantly with redshift, which will be investigated in a 
future work. A similar decline in met allicity has been 


established in the DLA pop ulation (e.g. Prochaska et al. 


2003a Rafelski et al. 2012) and has been interpeted as 


resulting from the ongoing enrichment of galactic ISM 
with cosmic time. Future work will perform a quantita¬ 
tive comparison between the two populations and explore 
the implications for the evolving enrichment of optically 
thick gas at z > 2. In passing, we emphasize the absence 
of low {XYH°} values at z « 2 which implies a reduced 
incidence of near-pristine gas with high H I columns at 
that epoch. 


6.3. Nucleosynthetic Patterns 

It is the conventional wisdom that LLSs primarily trace 
gas outside of the ISM of galaxies, e.g. within their dark 
matter halos (aka CGM) or at yet greater distances (iFu- 
magalli et al. 2011b Prochaska et al. 2013). Despite their 
separation from galaxies, we have demonstrated that the 
LLSs are generally enriched in heavy elements and pro¬ 
vided evidence that their metallicity frequently reaches 
~ 1/10 solar abundance. Therefore, a non-negligible 
fraction of this optically thick medium has been pro¬ 
cessed through the furnaces of a stellar interior and pre¬ 
sumably was transported from a galaxy via one more 
physical processes. One plausible transport process is 
an explosive event, e.g. a supernovae that expelled the 
gas shortly after enriching it. In this case, the gas may 
exhibit a distinct nucleosynthetic pattern from those ob- 


There is the possibility of a slight bias against our identifying 
metal-free SLLS but we have been as careful as possible to select 
systems based solely on the Lyo profile. 




Fig. 13. — Scatter plots of a series of ionic ratios, normalized to 
the solar relative abundances, against the A^hi values of the LLSs. 
The left two panels show ratios related to the o/Fe abundance. De¬ 
spite the preponderance of lower limits (especially for {0®/Fe+}), 
and concerns on the ionization corrections, we tentatively conclude 
that the LLSs exhibit super-solar a/Fe ratios, especially at large 
A^hi values. The measurements in the upper-right panel indicate 
that Si/C is possibly enhanced by a few 0.1 dex relative to solar 
although the majority of the sample is consistent with [Si/C] = 
0. Similarly, the {AR/C^} measurements are roughly consistent 
with the solar relative abundance or possibly sub-solar. In each 
panel, we indicate the expected offsets to the measurements that 
would be due to Type la (la) nucleosynthesis. Type II (II) nu¬ 
cleosynthesis, photoionization (PI), and differential depletion (D). 
Circles indicate a small or unknown impact. 


served for galactic ISM, i.e. if the supernovae ejecta did 
not mix prior to escaping the system. Additionally, the 
LLSs may couple the metal production within galaxies 
to the enrichment of the diffuse IGM (e.g. Aguirre et al. 
2001||Schaye et al.||2003||steidel et al.|)2U10l). 'This moti¬ 
vates comparison of the abundances tor these two diffuse 
and ionized phases. 

We may explore several ionic ratios that trace different 
nucleosynthesis channels. As with metallicity, one must 
account for i oni zation effects when interpreting the re¬ 
sults. Figure [T^ plots four pairs of ions from the dataset, 
again represented as {X®/Yl} with ionization corrections 
explicitly ignored. The figure also indicates the proba¬ 
ble offsets to the ratios if ionization effects were impor¬ 
tant, as estimated from photoionization calculations (e.g. 
Prochaska 1999). Similarly, we indicate the likely offsets 
from differential depletion and the dominant nucleosyn¬ 
thesis channels (Type la and Type II enrichment). 

The left-hand panels show two measures of a/Fe, a 
key diagnostic of the relative co ntributions of Type la 
and Type II SNe nucleosynthesis ( Tinsley||1979 ). Unfor¬ 
tunately, the {0°/Fe^} ratios are dominated by lower 
limits due to the saturation of O I 1302 and the non¬ 
detection of Fe II transitions. The values are nearly con¬ 
sistent with a solar abundance although there are at least 
two systems with {0°/Fe^} > -1-0.3dex suggesting an 
a-enhanced gas. Correcting for photoionization effects 
would only strengthen this conclusion. These two LLSs 
also exhibit a low metallicity ([0/H] « — 2) such that 
their chemical signature is very similar to that of metal- 
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poor Galactic stars ( McWilliain||1997 ). 

Turning to {Si~'’/Fe~'~}, the sample is dominated by 
measurements exceeding the solar abundance. This 
includes a non-negligible set of measurements with 
{Si^/Fe'*'} > +0.5 dex, and one may speculate that this 
represents the metal-enriched ejecta of Type II SNe. The 
{Si^/Fe’*'} ratio, however, is likely to require an ioniza¬ 
tion correction to accurately estimate Si/Fe. This could 
explain, in part, the positive {Si'''/Fe'''} values in Fig- 
On the other hand, the highest {Si^/Fe'*'} values 


ure 
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occur in LLSs with high A^hi values where one expects 
ionization effects to be minimap] We conclude, there¬ 
fore, that at least a subset of the KLS population exhibits 
super-solar a/Fe ratios indicative of Type II enrichment, 
even in higher metallicity gas. 

Previous studies of gas in the IGM at z ~ 2 have 


repor ted an enhanced Si/C abundance (Aguirre et al. 


20041. This result was derived statistically from the pixel 


optical depth method and is sensitive to the assumed 
mod el of the extr agalactic UV background (EUVB; see 
also Simcoej 20I1|). The results for LLSs offer a mixed 
picture (h'igure 13). There are a handful of positive 
{Si^/C+I values up to +0.5 dex, with the highest mea¬ 
surements at low metallicity. On the other hand, the 
sample is dominated by upper limits (from C II 1334 


saturation) and over half of these have {Si'''/C+} < 
+0.3 dex. Once again, photoionization corrections would 
only strengthen this result. As such, the LLS observa¬ 
tions do not appear to exhibit a high enrichment of Si/C 
than that previously inferred for the IGM. 


Figure [l3| also presents the set of {Al'''/C+} measure¬ 
ments that are not fully compromised by line-saturation. 
These data are consistent with the lighter element ratios 
in LLSs having solar relative abundances. The prepon¬ 
derance of upper limits, however, allows that Al could be 
under-abundant relative to C. 


6.4. Comparisons 

We have restricted the HD-LLS Sample to systems 
with A^hi < cm“^ to exclude the DLAs. This was 
partly motivated by the expectation that the majority of 
LLSs are predominantly ionized and therefore physically 
distinct from the neutral gas comprising DLAs. It was 
also motivated by the desire to examine this optically 
thick gas separately from the decades of research on the 
DLAs. Nonetheless, the A^hi = 10^°'^ cm“^ criterion is 
primarily an observationally defined boundary and one 
may gain insight into the nature of the LLSs through a 
combined comparison. Such analysis has been performed 
previously for the SLLS by “ilSom et al.| (|2013|). 

We consider two such comparisons here. Figure 
presents the {Si^/H*^} and {Fe~''/H°} measurements for 
the HD-LLS Sample t ogether with m e asurem ents from 
the sample of DLAs of|Rafelski et al. (20I2|). For both 
datasets, we have restricted to Zabs = [l.b, 3.3] to mini¬ 
mize trends related to redshift evolution. To zeroth or¬ 
der, the DLA measurements extend in a roughly continu- 


Such gas may also experience differential depletion, i.e. ele¬ 
vated Si/Fe ra tios in the gas phase to the refractory nature of these 
elements (e.g. |Jenkins||2009f. If the gas is predominantly ionized, 
however, the depletion levels may be modest and this effect would 
be small. 



log Vjn 


Fig. 14.— Comparison of the {Si+ZlF^} and {Fe+/H'^} mea¬ 
surements for the LLSs a nd DLAs. The latter are drawn from the 
abundance compilation of|Rafelski et al.|(|2012|l. For both ions, the 
DLAs show a continuous extension ol tne measurements observed 
in the LLSs. The right hand panels compare the distributions of 
the DLAs (darker) against those for the SLLSs (lighter), and em¬ 
phasize the commonality between the two datasets. In each case, 
we have treated upper and lower limit estimates as values. 


ous manner from the measurements of the LLSs. Indeed, 
comparing the samples of DLA measurements with the 
SLLSs (taking limits at their values), one observes over¬ 
lapping distributions with similar median values. The 
only notable distinction, perhaps, is the small set of LLSs 
with A^hi ~ 10^® cm“^ and high {XYH°} values (exceed¬ 
ing Odex for Si'*'). This suggests a higher incidence of 
highly enriched gas in the LLS, although we caution it 
could be partly an effect of ionization. The dispersion 
in the measurements is also larger for the LLSs, and is 
likely higher than suggested by the Figure given the pre¬ 
ponderance of upper/lower limits for the LLS/DLA. 

Turning to the higher ionization states. Figure 
presents the C’*'^ and Si"*"^ column densities from t 
LLSs and DLAs. Once again, the DLA distribution ex¬ 
tends in a nearly continuous manner from the upper end 
of the LLS data and the column density distributions for 
the SLLSs and DLAs are similar. Together, Figures 
and lend support to scenarios that envision LLSs as 
the outer layers of gas surrounding DLAs, i.e. these sys¬ 
tems frequently sample the same structures. Such physi¬ 
cal associations may be examined by studying DLAs and 
LLSs alo ng pairs of quasar sightlines, an active area of 
research (Ellison et al.|2007 Fumagalli et al.|20I4 Rubin 


et al.|[Ml^ . 


Examining the high-ion comparison further, there is at 
least one import distinction: the LLSs and especially the 
SLLSs show a much higher incidence of low AI(C+^) and 
A^(Si+^) values. This is unexpected given that (i) the 
LLSs trace highly ionized gas; (ii) the DLAs trace pre¬ 
dominantly n eutral gas that is physically distinct from 
the hig h-ions (Wolfe &: Prochaska||2000 Prochaska et al. 


2008a). The results presented here indicate that the gas 
layers giving rise to DLAs are embedded in a reservoir 
of highly ionized gas that frequently exceeds the typi¬ 
cal surface density in LLSs. This follows previous work 
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Fig. 15. — Comparison of the high-ion column densities mea- 
sured for the LLS and a representative se t of DLAs (drawn from 
|Rafelski et al.|2012| |Neelenian et al.|2013|l. Similar to the low-ion 
abundances, the measurements sliow a continuous transition from 
the LLS regime to higher Nm values. As such, there is substantial 
overlap in the distribution of measurements and limits (right-hand 
panels compare the values for DLAs [darker] and SLLSs [lighter]). 
One notable and surprising difference is that the SLLSs show a 
higher incidence of low columns of and Si"*"^. Despite trac¬ 
ing predominantly neutral gas, the DLAs also mark a reservoir of 
highly ionized gas that frequently exceeds the medium encompass¬ 
ing the LLSs. 


also correlated with A^hi- Examining ionic ratios sensi- 
live to the ionization state (e.g. C+^/C+, Si+^/Si+), we 
conclude that the LLSs are predominantly ionized with 
more highly ionized gas in lower Nm systems. 

Ratios of low-ion column densities to Nm indicate a 
wide spread in metal-enrichment within the LLSs, likely 
spanning four orders of magnitude. Only a small sub¬ 
set (< 5%) of the HD-LLS Sample have no positive de¬ 
tections of associated metals, consistent with primordial 
abundances. None of the LLSs with Nm > 10^®'^ cm“^ 
are ‘metal-free’. We conclude that a very high percent¬ 
age of high-density gas at z ~ 3 was previously enriched 
to > 1/1000 solar abundance. The HD-LLS Sample also 
exhibits a small subset 10%) of LLSs that have so¬ 
lar or super-solar enrichment. These likely represent the 
most enriched gas reservoirs in the high redshift universe. 

Lastly, we have examined several ionic ratios that are 
sensitive to the nucleosynthetic enrichment history of the 
gas. The preponderance of elevated Si+/Fe''" and 0°/Fe+ 
measurements suggest the LLSs have an a-enhancement 
characteristic of Type 11 nucleosynthesis. In contrast, 
the Si+ZC'*' and Al+/C“'' ratios are consistent with solar 
relative abundances. 

Future manuscripts on the HD-LLS Sample will: (i) 
study the metallicity distribution of the LLSs accounting 
for ionization effects and will estimate the contribution 
of optically thick gas to the cosmic metal budget; (ii) 
examine the kinematic characteristics to constrain the 
physical origin of the gas; (iii) offer constraints on the 
Nm frequency distribution for optically thick gas. 


which has inferred high q uantities of both neutral and 

m-ions 


Lehner et al. 
more closely 


ionize d gas for the DLAs (Fox et al 
20141. It further suggests that higK- 
trace higher density regions in the universe and/or may 
reflect a difference in the masses of the dark matter halos 
hosting LLSs and DLAs. 

Lastly, we have compared our measurements against 
the small set of literature values for SLLSs at z > 2 
(jDessauges-Zavadsky et al.||2003| |Som et al.||2013|) . The 
low-ion column densities are typically larger in those pub¬ 
lications, consistent with the higher Vhi values of the 
SLLSs that were sampled. 


7. SUMMARY 

We have constructed a sample of 157 LLSs at z ^ 2 —4 
observed at high-dispersion with spectrometers on the 
Keck and Magellan telescopes which constitute the HD- 
LLS Sample. In this manuscript, we present the complete 
sample and present column density measurements of H I 
and associated metal absorption. For the latter, analysis 
was restricted to transitions redward of the Lya forest 
and has focused on commonly detected species. These 
measurements and the associated spectra are made avail¬ 
able online with this publicatiorF] This constitutes, by 
roughly an order of magnitude, tEe largest high redshift 
sample of LLS analyzed in this manner. 

We have explored empirical trends in the column den¬ 
sity measurements and report statistically significant 
(> 99.99%) correlations between the low-ion (e.g. Si+, 
C+) columns and Ahi. High-ion species (Si+^, 0+^) are 
detected in nearly all LLSs and their column densities 
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APPENDIX 

APPENDIX: Measurements for Individual LLS 

Table lists measurements for all of the metal-line trans itio ns analyzed in this manuscript and Figures showing 
velocity ^ots are provided in the on-line materials (Figure [I^ shows one example). The analysis was restricted to 
lines outside the Lya forest and those lines that are not severely blended with another feature or compromised by 
sky-subtraction residuals. 
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TABLE 6 

IONIC COLUMN DENSITIES 


Quasar RA DEC z^bs ^HI A Ulim flg Nx Ion flg Njon <^i^) 

(J2000) (J2000) (A) (km/s) 


QOOOl-2340 

00:03:45 

-23:23:46.5 

2.18710 

19.65 

1334.5323 

-423,64 

2 

14.45 

99.99 

6,2 

2 

14.45 

99.99 






1335.7077 

-205,64 

4 

13.10 

99.99 










1548.1950 

-423,64 

0 

14.26 

0.01 

6,4 

1 

14.26 

0.05 






1550.7700 

-394,64 

0 

14.25 

0.02 










1302.1685 

-213,64 

0 

14.16 

0.04 

8,1 

1 

14.16 

0.05 






2852.9642 

-213,64 

4 

11.72 

99.99 

12,1 

3 

11.72 

99.99 






2796.3520 

-413,64 

2 

13.43 

99.99 

12,2 

1 

13.64 

0.05 






2803.5310 

-404, 64 

0 

13.65 

0.02 










1670.7874 

-405,64 

2 

13.00 

99.99 

13,2 

2 

13.00 

99.99 






1854.7164 

-213,64 

4 

12.40 

99.99 

13,3 

3 

12.40 

99.99 






1862.7895 

-213,64 

4 

12.69 

99.99 










1260.4221 

-399,64 

2 

13.81 

99.99 

14,2 

1 

13.75 

0.05 






1304.3702 

-423,64 

0 

13.55 

0.09 










1526.7066 

-399,64 

0 

13.83 

0.03 










1808.0130 

-213,64 

4 

14.86 

99.99 










1393.7550 

-411,64 

0 

13.78 

0.01 

14,4 

1 

13.74 

0.05 






1402.7700 

-421,64 

0 

13.64 

0.02 










1250.5840 

-79, 64 

0 

14.19 

0.13 

16,2 

1 

14.19 

0.13 






1608.4511 

-213,64 

4 

13.49 

99.99 

26,2 

1 

13.11 

0.05 






2344.2140 

-213,64 

0 

13.22 

0.07 










2374.4612 

-213,64 

4 

13.46 

99.99 










2382.7650 

-213,64 

0 

13.01 

0.04 










2586.6500 

-213,64 

4 

13.15 

99.99 










2600.1729 

-328,64 

0 

13.25 

0.04 










1317.2170 

-213,64 

4 

13.41 

99.99 

28,2 

3 

13.40 

99.99 






1370.1310 

-213,64 

4 

13.40 

99.99 










1454.8420 

-213,64 

4 

13.56 

99.99 










1741.5531 

-213,64 

4 

13.60 

99.99 










1751.9157 

-213,64 

4 

13.83 

99.99 










2026.1360 

-213,64 

4 

12.38 

99.99 

30,2 

3 

12.38 

99.99 

PX0034+16 

00:34:54.8 

-1-16:39:20 

3.75397 

20.05 

1548.1950 

-242,65 

0 

13.85 

0.02 

6,4 

1 

13.85 

0.05 






1550.7700 

-118,189 

2 

13.68 

99.99 










1670.7874 

-187,189 

0 

12.52 

0.04 

13,2 

1 

12.52 

0.05 






1854.7164 

-115,120 

4 

12.24 

99.99 

13,3 

3 

12.24 

99.99 






1862.7895 

-89,129 

4 

12.57 

99.99 










1526.7066 

-54,122 

2 

14.06 

99.99 

14,2 

2 

14.06 

99.99 






1808.0130 

-133,138 

4 

14.73 

99.99 










1393.7550 

-160,133 

0 

13.30 

0.02 

14,4 

1 

13.30 

0.05 






1741.5531 

-187,189 

4 

14.29 

99.99 

28,2 

3 

13.72 

99.99 






1751.9157 

-86,189 

4 

13.72 

99.99 










2026.1360 

-187,189 

4 

13.16 

99.99 

30,2 

3 

13.16 

99.99 


Note. — [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample.] 


Wolfe, A. M., Prochaska, J. X. 2000, ApJ, 545, 591 
Worseck, G., et al. 2014, MNRAS, 445, 1745 


Zafar, T., Popping, A., Sz Peroux, C. 2013, A&A, 556, A140 
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Fig. 16. — Velocity plots for the HD-LLS Sample 

Fig. Set 16. Velocity Plots for the Lyman Limit Systems 



























































































































